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ABSTRACT 

Tightest known quadruple systems VW LMi consists of contact eclipsing binary 
with Pi2 = 0.477551 days and detached binary with P34 — 7.93063 days revolving in 
rather tight, 355.0-days orbit. This paper presents new photometric and spectroscopic 
observations yielding 69 times of minima and 36 disentangled radial velocities for the 
component stars. All available radial velocities and minima times are combined to 
better characterize the orbits and to derive absolute parameters of components. The 
total mass of the quadruple system was estimated at 4.56 Mq. The detached, non- 
eclipsing binary with orbital period P = 7.93 days is found to show apsidal motion 
with J7 w 80 years. Precession period in this binary, caused by the gravitational 
perturbation of the contact binary, is estimated to be about 120 years. The wide mutual 
orbit and orbit of the non-eclipsing pair are found to be close to coplanarity, preventing 
any changes of the inclination angle of the non-eclipsing orbit and excluding occurrence 
of the second system of eclipses in future. Possibilities of astrometric solution and direct 
resolving of the wide, mutual orbit are discussed. Nearby star, HD95606, was found 
to form loose binary with quadruple system VW LMi. 
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1 INTRODUCTION 

The photometric variability of VW LMi (HIP 54003, 
HD 95660, sp. typ e F3-5V, y max=8.0), was found by the 
Hipparcos mission (|ESAlll997l L where it was correctly classi- 
fied as a W UMa-type eclipsing binary with an orbital period 
of 0.477547 days. The first ground-based photometric obser- 
vations of the system obtained in 1999 and 20 00 (taken in 
the B and V Johnson filters) were p ublished bv | Dumitresctil 
(|2000l ) . Analysis of these light curves (|Dumitrescrj2003l ) lead 
to the photometric mass ratio qph = 0.395 and inclination 
i = 72.4° and contact config uration for the system. Later 
iGomez-Forrellad et al.l (120031 ) presented new BV photom- 
etry and its preliminary analysis. Assuming convective en- 

* Based on the data obtained at the David Dunlap Observatory, 
University of Toronto 



velopes for both components and the temperature of the 
primary as Ti = 6700 K the authors estimated the mass 
ratio as g = 0.4, and inclination around 70°. Fourier analy- 
sis of it s Hip parcos light curve (hereafter LC) presented by 
iSelamI (|2004l ) yielded quite different parameters: q = 0.25, 
i = 72.5° and fill-out factor / = 0.4. The d iscovery of the 
second (non-eclipsing) binary in VW LMi bv lPribulla et all 
(2006) makes all previous photometric solutions almost use- 
less due to strong light contribution of the second pair of 
about (L3 -\-L4,)/{Li -\- L2) = 0.42 (at the maximum bright- 
ness of the contact pair) which was not taken into account. 

iPribuUa et al.l ([2006) presented long-term spectroscopy 
(209 spectra taken between 1998 and 2005) of the system ob- 
tained at David Dunlap Observatory (hereafter DDO) which 
enabled to disentangle all three orbits in this tight multiple 
system: the contact binary with the P12 = 0.4775 days pe- 
riod is orbiting another binary with P34 = 7.93 days in a 
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relatively tight, 355-days, mutual orbit. Using preliminary 
inclination angle of the contact-binary orbit found by photo- 
metric analysis, ii2 = 80.1°, the authors determined masses 
of all components and found that the orbits of the binaries 
are not coplanar. Light-time effect (hereafter LITE) with 
peak-to-peak range of 2A = 0.0074 days was predicted to 
be seen in the minima of the contact pair as a result of the 
mutual revolution of the binaries. The LITE wa s found in 
published minima by ISanchez-Baio et al.l l|2007h . The cor- 
responding orbital parameters, A = 0.0037(4) days, P1234 
= 353(2) days, e = 0.5(2) and uj — 2.8(4) rad, are rather 
preliminary due to few available minima. The eccentricity 
corresponding to their orbital solution is mu c h high er than 
predicted spectroscopically bv lPribulla et all (|2006i ). 

The spectral type o f VW LMi was estimated as F5V 
bv iPribulla et al.l (20061 ). while observed Tycho-2 {B - V) 
= 0.21 and 2MASS J-K = 0.34 colors correspond to F2V 
spectral type. Both determined spectral type and colors refer 
to the whole quadruple system. 

VW LMi is tightest quadruple system known 
I Tokovininll 20081 ) . Also it has the shortest period of the outer 
orbit within multiple systems harboring contact binaries. 
The ratio of the outer orbital period and orbital period of 
the non-eclipsing pair is only about P1234/-P34 ~ 44.5 hence 
we can expect secular orbital changes on the timescales 
as short as decades. The chances of resolving of compo- 
nents (binaries) by either speckle or long-baseline interfer- 
ometry are rather meager: the maximum angular separation 
of the components w as estimated to be only about lOmas 
IPribulla et al.ll2006h . Astrometric observations of VW LMi 
do not indicate its multiplicity. 

Goals of the present paper are as follows; (i) to present 
and analyze new photometric and spectroscopic observa- 
tions, (ii) to perform simultaneous solution of LITE using 
both minima times of the contact binary and radial veloc- 
ities (hereafter RVs) of the individual components, (iii) to 
assess possibility of the tidal disturbances of the inner orbits 
and resulting precession, (iv) determine absolute parameters 
of all four components. 



Table 1. Overview of telescopes and instruments/detectors used 
to obtain photometry of VW LMi. Abbreviations of the observa- 
tories (Obs.): Gl, G2 paviUons of the Stara Lesna Observatory, 
GSH - Grossschwabhausen observing station of the Jena Univer- 
sity, KO - Astronomical Observatory at Kolonica saddle (belongs 
to the Vihorlat Observatory), RO - Roztoky Observatory. Instru- 
ments using a photolense at KO and RO are named using its focal 
length. 



Obs. 


Telescope 


Detector 


Filters 


Gl 


508/2500 


SBIG STIO-MXE 


UBV{RI)c 


G2 


600/7500 


EMI 9789QB 


UBVN 


GSH 


250/2250 


CTK102 


B 


KO180 


60/180 


MEADE DSI Pro 


N 


KO400 


80/400 


MEADE DSI Pro 


VN 


KO 


300/2400 


SBIG ST9-XE 


B 


RO200 


100/200 


MEADE DSI Pro 


N 



Table 2. Parameters of the variable and comparison stars used. 
The proper motions and (B —V) = 0.8 5(-B— V)r colors were taken 
from the Tycho-2 Cat alog |H0g et al . 2000), parallaxes from the 
Hipparcos C atalogue iIESAI 19971). infrared colors, (J — K), from 
the 2MASS l|Skrutskie et al.ll2006h . Radial (RV) and rotational 
velocities (iisini), and spectral types were determined from the 
DDO spectroscopy; for VW LMi the listed RV corresponds to the 
systemic velocity of the whole quadruple. 





VW LMi 


cmpl 


cmp3 


GSC 


2519-2347 


2519-1195 


2519-1326 


HD 


95660 


95606 


95527 


HIP 


54003 


53969 




/XqCos<5 [mas.y~^] 


12.7(11) 


10.7(11) 


-9.8(12) 


^1S [mas.y-l] 


-5.0(12) 


-3.4(11) 


-67.9(12) 


RV [km.s-i] 


-0.15 


7.80 


11.10 


V sin i [km.s^-*-] 




32 


< 13 


TT [mas] 


8.04(0.90) 


7.07(1.25) 




[B-V) 


0.340(21) 


0.398(28) 


0.570(34) 


(J-K) 


0.208(30) 


0.266(31) 


0.300(35) 


sp. type 


F3-5V 


F5V 


GOV 



2 NEW OBSERVATIONS 

2.1 Photoelectric and CCD photometry 

To secure as many minima times as possible, photometry 
was performed at several observatories in Slovakia and Ger- 
many. Since the system was too bright for most of the in- 
struments equipped with CCD cameras, photolenses were 
extensively used and on larger telescopes, the observations 
were mostly performed in the Johnson B filter which enabled 
reasonable exposure times. The information on individual 
observatories and instruments used is given in Table [T] pa- 
rameters of variable and comparison stars are comprehen- 
sively listed in Table [2] 

Photoelectric photometry of VW LMi was obtained us- 
ing 60cm Cassegrain telescope in the G2 pavilion of the 
Astronomical Institute of the Slovak Academy of Sciences. 
VW LMi was observed in the Johnson BV filters and with 
neutral filter (optical glass BK7, denoted in tables as A''). 
BV observations were transformed to the international pho- 
tometric system. The magnitude differences with respect to 



the comparison star HD95606 were corrected for the differ- 
ential extinction using nightly or seasonal extinction coef- 
ficients. Since the angular distance of the comparison and 
variable star is only 5.65 arcmin, the differential extinction 
correction was rather negligible, and never exceeded 0.002 
mag. Occasionally observed check star, HD95527, showed 
good stability of the comparison star, HD95606. 

Most of the CCD photometry was secured using 50cm 
Newtonian telescope in the Gl pavilion of the Astronom- 
ical Institute of the Slovak Academy of Sciences. Practi- 
cally all observations were taken through the Johnson B 
filter, only on March 25, 2007 UBV{RI)c photometry was 
performed. The data were left in the instrumental system 
since the major purpose was securing minima instants. Fur- 
ther large set of data was obtained using various small 
telescopes/photolenses at the Astronomical Observatory at 
Kolonica saddle (48° 56' 6" N, 22° 16' 26" E). Several times 
of minima were also obtained at the Roztoky Observatory 
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HJD 2 453 825.72206 
0.770 



HJD 2 454 180.67292 
0.044 




HJD 2 454 605.58769 
0.824 



-400 -200 200 400 -400 -200 200 400 -400 -200 200 400 

RV (km/s] RV (km/s) RV(km/s) 

Figure 1. Broadening functions of VW LMi. From left to right: 
all components visible, contact binary close to the conjunction, 
and detached binary close to the conjunction. Original broadening 
functions were smoothed convolving with a Gaussian function (cr 
= 10 km s~^). Heliocentric Julian dates and phases of the contact 
binary are given in individual panels. 



(49°33' 57" N, 21°28' 54"E). CCD photometry of VW LMi 
was also obtained at Grossschwabhause n observing statio n 
of the Jena University Observatory (see iMugrauerl l|2008l )). 
CCD reduction was carried out by usual way (bias, dark, 
flat-field corrections) and aperture photometry was per- 
formed by procedures in MIDAS environment (only at Stara 
Lesna Observatory), and C-Munipack packag^B- HD95606 
served as comparison star at all observatories. 



2.2 Medium dispersion spectroscopy 

New spectroscopic observations were obtained using the slit 
spectrograph in the Cassegrain focus of 1.88m telescope 
of DDO. The observations were performed between March 
2007 and very beginning of July 2008, when DDO ceased 
to operate. The spectra were taken in a window of about 
240 A around the Mg I triplet (5167, 5173 and 5184 A) 
with an effective resolving power of J? = 12,000 - 14,000. 
The diffraction grating with 2160 lines/mm was used with 
the sampling of 0.117 A/pixel. One-dimensional spectra were 
extracted by the usual procedures within the IRAF environ- 
men10 after the bias subtraction and the flat-field di vision. 
Cosm ic-ray trails were removed using a program of IPvchI 
(|2004l ). 

All exposures were 900 seconds long and lead to S/N 
of about 100-150. The observations were analy zed using 
the te c hniqu e of broadening-functions, described in lRucinskil 
j 19921 . 12OO2I ). For VW LMi we used HD222368 as the tem- 
plate (F7V, vsini — 3 km s~^). This star provided the best 
match to the observe d spectra - integral of BF was always 
0.95 - 1.0(0. Similar to lPribuUa et al.1 (I2OO6D we firstly fitted 
simple multi-component Gaussian model to the extracted 



^ http://munipack.astronomy.cz/ 

^ IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the NSF. 

^ Late spectral type of the best template indicates later spectral 
type of the second binary in VW LMi 



BFs (see Fig. [TJ, then detached pair with slowly-rotating 
components was subtracted. The wide features correspond- 
ing to the contact pair were subsequently modelled by the 
rotational profiles which effectively corresponds to approxi- 
mation of the components of the contact binary by rotating 
limb-darkened spheres. Since shapes of the profiles slightly 
depend on the li mb darkenin g we fi xed u = 0.53 according 
to the tables of Ivan Hammd l| 19931 ) appropriate for main- 
sequence F3V star at 5184 A. RVs determined from the new 
spectra are given in Table |3] Several BFs extracted from ob- 
servations taken close to conjunction of the contact binary 
did not allow to determine their RVs. In case of non-eclipsing 
binary, RVs of its components could be reliably determined 
from all new spectra. 



3 MINIMA DETERMINATION 

Most widespread approach to obtain instants of minima of 
eclipsing binaries is to use K wee & van Woerden's method 
(|Kwee fc van Woerden|[l95(j ). From our experience the er- 
rors estimated using their formula (14) are often unreal- 
istically small. The real uncertainties are very probably 
dominated by systematic errors. The principal problem in 
case of CCD photometry is scattered light which cannot 
be fully corrected by flatfielding. This results in spurious 
shifts/trends in differential photometry in case that the field 
is not perfectly guided making times of minima systemati- 
cally shiftecQ. Our photometry could possibly be improved 
by using an algori thm based on Princip al Component Anal- 
ysis proposed by iTamuz et al.l (|2005l ). Unfortunately, the 
frames were obtained at several observatories with different 
setups and even different orientation of the field. Systematic 
errors in minima positions were partially removed by fitting 
technique proposed below. 

Since LC of the system appears to be very stable (we do 
not see any asymmetries or changes) fitting templates were 
prepared to obtain instant of conjunction (minimum) for any 
sufficiently long photometric sequence. Such a way we made 
use not only of the minima but of other LC segments. The 
template LCs were represented by symmetric trigonometric 
series of the 10th order. Even if the eclipsing pair is a con- 
tact binary, amplitude of its LC depends on the wavelength. 
Using CCD photometry of March 25, 2007 we see the follow- 
ing amplitudes AB = 0.46, AV = 0.43, ARc = 0.41, and 
Ale = 0.38. The differences in amplitude primarily result 
from wavelength-dependent limb darkening and hght con- 
tribution of the third component (see Section [5]). Due to the 
differences in filter transparencies and wavelength response 
of detectors we had to form a template LC for each filter 
separately and the fitting LC was scaled to match the obser- 
vations. We also noted small (~ 0.02 mag) shifts of the LCs 
observed even with the same instrument. Sometimes even 
slight slopes of the LC were recorded. These shifts/slopes 
are very probably caused by scattered hght combined with 
different pointing of the telescope. 

To obtain good fits of the observed LCs by templates 
T{x) we formed the following fitting function: 



* Unfortunately, neither of the telescopes uses an autoguider for 
pointing 
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Table 4. New times of minima of VW LMi determined by the template fitting to observed LCs. Observations in individual filters were 
treated separately. If both minima were observed in given night (indicated in column "Type" ) the instant given in the table refers to the 
first one. The standard errors of the last digit are given in parentheses. 



HJD 


Fil. 


Obs. 


Type 


HJD 


Fil. 


Obs. 


Type 


HJD 


Fil. 


Obs. 


Type 


2 400 000+ 








2 400 000+ 








2 400 000+ 








53461.45927(21) 


B 


G2 


I 


54063.65646(14) 


B 


Gl 


I 


54245.35634(11) 


B 


Gl 


II 


53461.45883(16) 


V 


G2 


I 


54066.52319(13) 


B 


Gl 


I 


54415.61257(30) 


V 


KO 


I 


53461.45969(20) 


N 


G2 


I 


54085.62425(21) 


V 


KO400 


I 


54469.57672(29) 


B 


Gl 


I 


53465.51832(22) 


B 


G2 


II 


54088.49087(24) 


B 


Gl 


I 


54469.57628(29) 


V 


Gl 


I 


53465.51792(15) 


V 


G2 


II 


54095.65269(10) 


B 


Gl 


I 


54469.57549(20) 


N 


KO180 


I 


53465.51843(22) 


N 


G2 


II 


54114.51541(08) 


B 


Gl 


II 


54499.66040(07) 


B 


GSH 


I 


53767.57388(15) 


N 


G2 


I 


54116.66409(07) 


B 


Gl 


I 


54505.38874(26) 


N 


KO400 


I 


53791.44941(26) 


B 


G2 


I 


54117.61952(14) 


N 


KO400 


I 


54509.44938(43) 


N 


KO400 


II 


53791.44924(21) 


V 


G2 


I 


54148.65937(14) 


B 


Gl 


I 


54512.55203(10) 


B 


GSH 


I 


53791.45004(33) 


N 


G2 


I 


54148.65953(35) 


V 


KO400 


I 


54521.62540(08) 


V 


GSH 


I 


53794.31433(08) 


B 


Gl 


I 


54149.61531(13) 


B 


Gl 


I 


54532.36983(15) 


N 


KO400 


n,i 


53797.41955(16) 


B 


Gl 


II 


54162.50371(27) 


V 


KO400 


I 


54556.48551(26) 


V 


KO400 


I 


53802.43347(31) 


B 


G2 




54167.28101(36) 


V 


KO400 


I 


54557.43985(19) 


N 


RO 


I 


53802.43276(22) 


V 


G2 




54172.53421(10) 


B 


Gl 


I 


54571.28992(25) 


N 


RO 


1,11 


53802.43260(24) 


N 


G2 




54173.48896(12) 


B 


Gl 


I 


54581.31655(14) 


N 


RO 


i,n 


53814.37177(07) 


B 


Gl 




54173.48891(17) 


V 


KO400 


I 


54582.51234(24) 


N 


RO 


II 


53815.32588(05) 


B 


Gl 




54175.39958(16) 


V 


KO400 


I 


54584.41972(65) 


N 


KO400 


II 


53815.32524(08) 


N 


G2 




54185.42614(26) 


B 


Gl 


I 


54584.42043(18) 


N 


RO 


II 


53833.47072(07) 


B 


Gl 




54185.42681(32) 


V 


Gl 


I 


54585.37523(15) 


N 


RO 


II 


53845.40960(09) 


B 


Gl 




54185.42691(41) 


R 


Gl 


I 


54586.32887(24) 


B 


GSH 


II 


53850.42278(11) 


B 


Gl 


II 


54185.42709(39) 


I 


Gl 


I 


54588.47856(17) 


B 


GSH 


I 


53867.37715(08) 


B 


Gl 




54185.42676(24) 


V 


KO400 


I 


54593.49393(18) 


N 


RO 


II 


53894.35622(16) 


B 


Gl 




54189.48500(30) 


V 


KO400 


II 


54594.44843(18) 


N 


RO 


II 


54026.64630(10) 


B 


Gl 




54191.39733(21) 


V 


KO400 


II 


54595.40379(19) 


V 


KO400 


II 


54027.60191(19) 


B 


Gl 




54195.45500(15) 


V 


KO400 


I 


54599.46132(19) 


B 


GSH 


I 


54057.68914(14) 


B 


Gl 




54212.40823(10) 


B 


Gl 


II 


54615.46100(25) 


N 


RO 


II 


54059.59910(10) 


B 


Gl 




54213.36179(07) 


B 


Gl 


II 


54616.41716(20) 


N 


RO 


II 



F{x) =A + B.x + C.T{x ~ D). 



(1) 



2006 



2007 



2008 



This allowed for shifting, scaling and "slanting" of the 
template LC. Fixing of the parameters was judged according 
to the appearance of individual LCs. All new times of min- 
ima are listed in Table |4] The formal errors given are still 
much smaller than systematic errors. Times of "minima" de- 
termined from sections of LC not containing minima were 
omitted from further analysis because of significantly higher 
errors. 

The resulting (0-C) diagram including published min- 
ima times is shown in Fig. (2] The minima times of VW LMi 
can b e found in [Go mcz-ForrcUad ct al. ( 2003jj_ Dumitrescul 
20031) : iDerman fc Kal ci (2003): [Porowskj liooBI) : iDworakI 
20051 . l200d ): iHubsTher, (.2007. ): lNeiso3~^Ooland there 
are four unpublished CCD times of minima observed by 
K. NakajimsQ - HJD 2 453 461.9371(2), 2 453 462.1758(2), 
2 453 886.0001(3), and 2 454 603.9984(2). The residuals, from 
rather arbitrary ephemeris given in the figure, clearly show 
wave-like behavior with one-year periodicity. 




-0.010 



53400 53600 



53800 54000 54200 
JD - 2400000 



54400 54600 



Figure 2. (0-C) diagram of VW LMi showing best covered in- 
terval since 2005. The (O-C) data were given with respect to the 
following ephemeris for the primary minimum: 2 452 500.1510 + 
0.477552 xE. The primary and secondary minima were plotted us- 
ing the same symbols. The optimum fit corresponding to Table |5] 
is plotted. 



4 GLOBAL FIT OF THE DATA 

Multiple-dataset fitting is rather widespread technique 
nowadays. It is typical to combine RVs and visual orbits 



see |http://www. kusastro.kyoto-u.ac.jp/vsnet/ 1 



or RVs and inter fe rometric visib il ities o f binary stars (e.g., 
iPourbaixl (|l999t ). iTango et al.] (|2006|)). Eclipsing-b inary 
modelling programs (see Ivan Hamme fc WilsonI (|20o3)) en- 
able to combine RV data LCs to determine third-body pa- 
rameters. The case of VW LMi is, however, special and un- 
precedented, requiring ad-hoc approach. 
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Table 3. New radial velocities for all four components of the 
quadruple system VW LMi. While radial velocities of the de- 
tached pair (RV3,RV4; P = 7.93 days) were determined by multi- 
component Gaussian fitting to the observed BFs, radial velocities 
of the contact binary (RV1,RV2; P = 0.4775 days) were deter- 
mined by rotational-profile fitting to the residual BFs after re- 
moval of the detached pair from BFs. Radial velocities RVl, RV2 
are not given for observations performed close to the conjunc- 
tions. Heliocentric Julian dates were not corrected for the LITE 
on the 355-days orbit. Components denoted as "1" and "3" are 
more massive ones in either of the binaries. 



HJD 


RVl 


RV2 


RV3 


RV4 


2 400 000-1- 


[km/s] 


[km/s] 


[km/s] 


[km/s] 


53823.67103 






-16.88 


62.46 


53825.72206 


84.55 


-245.57 


68.77 


-29.03 


53825.73278 


78.31 


-237.41 


69.77 


-28.65 


53825.74469 


78.44 


-230.85 


69.92 


-29.05 


53825.75545 


72.69 


-206.33 


70.55 


-28.90 


53835.80611 






53.04 


-15.75 


53836.77953 






3.88 


36.12 


53836.78848 






3.86 


36.66 


54167.69771 


68.79 


-188.87 


87.41 


-39.54 


54167.70691 


58.05 


-180.22 


87.23 


-39.55 


54177.69983 


71.34 


-231.89 


12.57 


34.40 


54179.77080 


-112.25 


173.30 


-40.74 


84.99 


54180.67292 






-11.94 


57.31 


54182.66132 


— 114.17 


203.01 


—27.91 


73.79 


54183.65132 


-118.18 


213.28 


-42.67 


83.14 


54189.64139 


72.54 


-210.75 


34.34 


6.43 


54190.66416 






-33.23 


71.64 


54193.67997 


-108.38 


231.24 


37.19 


3.20 


54193.75446 






-0.17 


41.89 


54200.62450 


81.57 


-220.56 


49.60 


-17.90 


54507.71485 


58.29 


-204.99 


69.31 


-26.19 


54515.82080 


64.89 


-205.72 


75.41 


-31.01 


54518.69236 


61.91 


-195.11 


13.48 


31.94 


54521.78963 


-99.01 


170.49 


-4.12 


56.24 


54542.78635 


-104.76 


193.21 


-6.60 


44.71 


54546.74822 


70.95 


-166.15 


43.40 


-5.20 


54550.69995 


67.16 


-184.81 


-4.91 


43.48 


54605.58769 


101.56 


-217.66 


7.86 


-10.86 


54619.58209 


-67.07 


187.80 


56.99 


-68.57 


54626.64631 






38.25 


-55.24 


54626.65701 






38.87 


-55.73 


54628.60047 






34.47 


-55.86 


54638.59131 






-65.05 


40.82 


54642.60165 


-79.75 


234.95 


37.10 


-62.29 


54647.59897 


121.47 


-208.31 


-77.27 


51.63 


54649.58784 






-2.21 


-28.02 



Since our photometry was performed at several obser- 
vatories, often without filter, we decided not to fit a global 
model to LCs together with RVs and to use only minima 
times to better characterize the outer orbit. Rigorous mod- 
elling of those datasets (four RV curves and times of min- 
ima) is quite complex since times of all RV observations 
should be properly corrected for the LITE caused by the 
mutual revolution of the binaries. This effect is minor in 
the non-eclipsing pair, where amplitude of LITE is about 
0.004 days, which amounts to only about ±0.0005 in phase. 
In case of eclipsing pair the amplitude of LITE is about 
0.0037 days, and P12 = 0.4775 days make for ±0.0077 phase 



shifts which cannot fully be neglected. In fact, correcting ob- 
served times of RV measurements for the LITE significantly 
decreased global ■ LITE is also manifested directly in ob- 
served times of minima. The informational contents of the 
minima is, unfortunately, very similar to having observed 
only RVs of the eclipsing pair. Integrals of RV curves are 
radial distances from the observer corresponding directly to 
the LITE through the finite speed of light. Both datasets are, 
however, nicely complementary: during the spectroscopic 
conjunction, when the RVs of the components are not mea- 
surable, we can observe minima giving the radial distance of 
the eclipsing binary in the outer orbit. Moreover, the precise 
minima times are easy to observe using a small telescope or 
photolense equipped with a cheap CCD camera. As we will 
see later, the relative orbit of mass center of the contact bi- 
nary is much better defined by the LITE than RVs. This is 
crucial to establish mass ratio of binaries, inclinations and 
masses of components. 

In further description we will u se similar notation of 
parameters as in lPribuUa et al.l l|2006l ): common parameters 
(orbital period, eccentricity, time of periastron passage, lon- 
gitude of periastron) for orbits will be denoted by index "12" 
(contact binary), "34" (detached non-eclipsing binary) and 
"1234" (mutual orbit). Semi-amplitudes of RV changes for 
the individual stars will be denoted by separate index "1", 
"2", "3", or "4" while semi-amplitudes of systemic velocity 
changes of both binaries as whole in the mutual orbit will 
be denoted as "12" and "34". Finally, systemic velocity of 
the mass center of the whole quadruple as Vq. For RVs of 
components of contact binary we can write: 

RVi = Vo ± A'l2[ei234 COSiJl234 ± COs(!yi234 ± ^^1234)] ± 

+ (-l)'"'"^A'4ei2 COSCJ12 ±cos(!/i2 +0.^12)] (2) 

where index i = 1, 2 denotes the component. For the 
detached binary we similarly have: 

RVi = Vo - A'34[ei234 COS(i)l234 ± COs(i/l234 ± ^^1234)] ± 

+ { — iy^^ Ki[e34 COSUJ34. + COs(!/34 ± t^34)]. (3) 

with i = 3, 4. True anomalies of inner orbits {1^12, va) 
have to be corrected for the LITE caused by the mutual 
revolution of binaries. If we accept the plane parallel to the 
sky and intersecting mass center of the quadruple system as 
the reference plane, then LITE seen in either of the binaries 
is: 

^rp _ -K"j-Pl234(l — £1234)^'^^ Sin(i/i234 ± ^1234) 
^ 27rc 1 + 61234 cos i^l234 ' 

where index j is either "12" or "34". In case of the 
eclipsing pair, where the orbital eccentricity is zero, the pre- 
dicted times of minima occurring at epoch E would be: 

Trmn = Tl2 ± ^12^ ± QE'^ ± AT12 (5) 

if true anomalies in the eclipsing binary orbit are 
counted from the upper conjunction of the more massive 
component (then a;i2 = ■k/2 with ei2 = 0). 

While amplitude of the RV changes of the systemic ve- 
locity of the detached-pair orbit around the common center 
of gravity is rather reliable due to very well defined sharp 
peaks in BFs [v sin i < 12 km s"'^ which is limit given by the 
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spectral resolution), the RVs of contact-binary components 
are rather imprecise due to a "cross-talk" with the second 
binary and wide profiles. This makes size of the relative or- 
bit of the contact binary (related to K12) rather unreliable. 
The shape of the orbit is, fortunately, coupled through well 
defined orbit of the detached pair. The size of the contact- 
binary orbit is much-better defined by LITE seen in the 
observed minima. 

The best fits to the data (4 RV curves and the LITE 
orbit) were performed using differential corrections method 
using analytic derivatives of the functions. The dependence 
of true anomaly on mean anomaly (M) and eccentricity was 
represented in derivatives by the truncated series up to the 
second degree in orbital eccentricity: 



5 9 

^ M + 2e sin M + -e sin 2M. 

4 



(6) 



Since all orbits are being close to circular this approx- 
imation was sufficient to make the optimization process 
quickly convergent. The sum of normalized was used as 
the merit function: 



S = 



E 

.7 = 1 



5 

J = l 



(7) 



where j is index of the dataset, Uj is number of ob- 
servation for the given dataset, i is index of individual ob- 
servation and aj is average standard deviation of individ- 
ual datapoint for the j-th dataset. Average uncertainty of a 
datapoint for individual datasets was estimated fitting sep- 
arately times of minima and RVs. For minima times we got 
o'min = 0.0008 day (which is significantly more than typical 
error in Table |4)|, for RVs anvi = 11.3 km s~^, aiiv2 = 
13.1 km s~^, aRV3 ~ 2.6 km s~^, and aRV4 = 2.2 km s~^. 
The estimated average uncertainties were found to be realis- 
tic since normalized for individual datasets were close to 
unity. Observations within individual datasets were assigned 
the same errors (as given above) due to the fact that both 
times of minima and RVs are mostly affected by system- 
atic errors which cannot reliably be determine d. One mini- 
mum (2450,877.4658), published bv iGomez-Forrellad et al.l 
l|2003h . giving very large (O-C) (probably due to a typing 
error), was omitted from further analysis. 

Fitting theoretical curves to five datasets required 18 
parameters. It is interesting to note, that LITE of the eclips- 
ing pair required only one additional parameter - Q, due to 
the fact that the orbital period of the contact binary was 
found to be continuously increasing. Since orbital eccentric- 
ity of the detached binary was found to be only 634 = 0.03, 
we also considered a circular orbit. This, however, lead to 
substantially worse fit to the data. Resulting parameters are 
given in Table [S] Corresponding fits are shown in Figs. [5] and 

El 

Unfortunately the orbital period of the mutual orbit is 
close to one year. The system can be observed for minima 
from beginning of November till beginning of June. The ob- 
serving interval to get spectrum just before dawn and just 
after twilight is about a month longer. Therefore, to cover 
the whole outer orbit, the system would have to be moni- 
tored at least 12-15 years (for ground-based observations). 
Until then parameters depending on overall shape of the 
RV curve (all except orbital period) cannot be reliably de- 
termined. This affects reliable determination of mass ratio of 



Table 5. Simultaneous fit to radial velocities of both binaries 
and times of minima of the eclipsing pair. The designation of 
parameters is as follows: the index "12" refers to the orbit of 
the contact pair, while the index "34" refers to the orbit of the 
second, detached binary. Parameters of the mutual orbit of these 
binaries are indexed as "1234" . The orbit of the contact pair is 
assumed to be circular (ei2=0, LU12 = Tr/2). Masses of the systems 
are abbreviated as M1234 = Mi + M2 + M3 + M4, M12 = Afi -|- 
M2, M34 = M3 + M4,. Standard errors of parameters are given in 
parentheses 

Contact (eclipsing) pair - circular orbit 



P12 [days] 
Q [days] 
T12 [HJD] 
Ki [km s-i 
K2 [km s-i 
M12 sin^ ii2 

x'^CRVi) 

X2(RV2) 
2(MIN) 



[Me 



X 



0.47755106(3) 
1.63(9) 10-10 
2,452,500.1497(2) 
105.8(1.0) 
250.2(1.2) 
2.231(23) 
1.086 
1.058 
1.086 



Detached (non-eclipsing) pair 



P34 [days] 
634 

0)34 [rad] 
T34 [HJD] 
K3 [km s-i] 
K4, [km s-i] 
M34sin3i34 [M 
X2(RV3) 
x2(RV4) 



©J 



7.93063(3) 
0.035(3) 
1.90(9) 
2,452,274.54(11) 
63.99(23) 
65.53(27) 
1.785(11) 
0.862 
0.851 



Mutual wide orbit 




-P1234 [days] 


355.02(17) 


61234 


0.097(11) 


t^l234 [rad] 


2.20(12) 


T1234 [HJD] 


2,453,046(6) 


ftTia [km s"!] 


21.61(49) 


K34 [km s-l] 


23.22(33) 


M1234 sin^ ii234 [Mq] 


3.32(10) 


Vb [km s-i] 


-0.15(25) 



the binaries M34 /M12 = K12 / K34 and projected total mass 



(Ml + M2 + M3+ M4) sin'' 11234. 



5 LIGHT-CURVE SOLUTION 

The eclipsing pair revolves in rather tight orbit around an- 
other binary which results in phase shifts of about ±0.008 
which cannot be neglected as seen in long-term modelling 
of RVs (see Section U]). Therefore, we used photometry from 
a short interval compared to mutual, 355-days orbit. To de- 
termine photometric elements of the eclipsing pair we used 
BV photoelectric LCs obtained on March 31 and April 4, 
2005 at the Stara Lesna observatory giving full coverage of 
the phases in relatively short time compared to the mu- 
tual orbital period (about 1.2%). The orbital period was 
fixed at 0.477551 days but time of upper conjunction of 
the more massive component was adjusted. The third light, 
(L3 -I- L4,)/{L\ + L2) = 0.42, is known from spectroscopy. It 
was derived from a 240A section of the spectrum centered 
at 5184A fairly close to the Johnson V filter. General expe- 
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Figure 3. Radial velocities of all four components and their best 
fits assuming hierarchic quadruple model. Phases of the measure- 
ments in inner orbits were corrected for LITE in the outer orbit. 
Longitude of the periastron of the detached pair is assumed to be 
constant. Top panel and center panels show RVs of the contact 
and detached pairs, respectively, after correction for the revolu- 
tion in the outer 355-days orbit. The bottom panel displays veloc- 
ities of individual components corrected for motion in inner orbits 
- reflecting mass-center revolution (observations of the secondary 
component of the contact binary were not plotted for clarity). 
Because the outer orbit period is close to one year, the seasonal 
observing interval only slowly moved in the orbital phases be- 
tween 1998 and 2008 and all phases have not been covered yet 
(bottom). The best fits corresponding to the global solution (Ta- 
ble [Sjl are plotted. 
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Figure 4. Best fits to the BV LCs for the adopted solution with 
i = 79.0° fixed. 



rience with BF techn ique is, however, that the t hird hght is 
often overestimated (|Rucinski fc Pribulla|[2008l ) due to the 
continuum normalization problems of contact binaries and 
difference of spectral types of additional components. 

Analysis of LCs was performed using code ROCHE de- 
scribed in .Pribulla (2004). Since eclipsing pair is a contact 
binary, circular orbit and synchronous rotation (Fi = F2 = 
1) was assumed in al l fits. The model-atmosphere grids (see 
iLeieune et al.l (|l997l )') were selected for the solar metallic- 
ity. The temperature of the primary component was fixed 
according to its spectral type, F5V (P ri bulla o t al. 20061), 
at 6460 K using calibration of IPopperi l|l980l ) . Convective 
transfer of energy in the common envelope was assumed 
{Al = A2 = 0.5,(71 = (72 = 0.32), and limb-darkening 
was automatically interpolated from the extensive tables of 
[van Hammg (1993,) for mean logg and T^ff of the compo- 
nents. Due to the fact that exact values of third hght for 
B and V passbands are unknown, unlike in usual approach, 
we performed several solutions with the orbital inclination 
fixed between 75 to 85°. The mass ratio was not adjusted 
since it tightly correlates with orbital inclination, and third 
light, but was rather fixed at spectroscopically determined 
value Qsp = 0.423. Optimization of the parameters lead, as 
expected, to contact configuration for every trial inclina- 
tion. Surprisingly, third hght was found to be almost iden- 
tical in B and V passbands, indicating very similar spectral 
type of the second non-eclipsing pair. For i — 79.0°, third 
light in both passbands was close to spectroscopically de- 
termined value L34/Li2{B) = 0.41±0.03 and L34/ii2(V) = 
0.42±0.03. The resulting fill-out factor is rather hight, F = 
0.63±0.05. The best fits to the data are shown in Fig. H 
The major deficiency of the fits is insufficient interpretation 
of the primary minimum, which is sharper than predicted 
suggesting either higher mass ratio, surface inhomogeneities 
or deficiency of the Roche model for this system. It cannot be 
excluded that the observed LC is affected by some night-to- 
night changes resulting in systematically wrong parameters 
and fits. 
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6 LONG-TERM EVOLUTION OF ORBITS 

The ratio of the orbital periods of the outer orbit and orbit of 
the non echpsing pair is P1234/-P34 = 44.76. Corresponding 
period ratio in the case of the contact binary is as large as 
P\234,/Pi2 = 743. Therefore, visible precession of the orbit 
can be detected (on human timescales) only in the case of 

the detached pair. ^ 

According to ISoderhielml (|l975l ) the nodal precession 
period caused by the perturbations in the close triple is: 

_ 4 /Mi234 \ ^"1234/, 2 n3/2 il234 

rnode. — - I —77 — f5 U — 61234 J 77 \ — r \ 

3 V M34 / P34 (1/34 + il234) COSJ 

where Z/34 is angular momentum in the detached binary 
orbit, and Z/1234 is (vector) sum of angular momenta in the 
detached binary orbit and the outer orbit. In our case we 
can neglect the term in the last parentheses since (I/34 + 
ii234)/-£'i234 ~ 1.08 and orbits are close to being coplanar; 
hence cos j ~ 1.00. Then we get Pnode ~= 120 years. Since 
the detached pair on P34 — 7.93-days orbit is non-eclipsing 
binary (at least at present) we can detect possible precession 
of the orbit only from the RVs. Precession of the orbital 
plane should results in changes of the amplitude of RVs and 
apsidal motion. 

Possibility of long-term changes of the 7.93-days orbit 
was assessed by dividing RV observations into individual 
observing seasons. Since LITE in the eclipsing pair is not 
affected nor RVs of its components, simultaneous fits in sep- 
arate time intervals used all measurements for the eclipsing 
pair with the parameters of the contact-binary orbit and 
outer orbit fixed to those of the global solution. 

Five solutions given in Table [6] can be de- 
scribed/interpreted as follows: (i) eccentricity is stable - so- 
lutions are reliable and small eccentricity is not spurious (ii) 
semi-amplitudes of RVs remain stable within the errors - 
mutual orbit and orbit of the detached pair are very prob- 
ably close to coplanar (iii) there is definite apsidal motion. 
The coplanarity of the orbits means that the orbit of the 
non-eclipsing binary very probably not changes its inclina- 
tion angle and we cannot ex pect ap p earan ce of the eclipses 
in this binary in future (see iMaven (|2005| )). The period of 
the line of apsides motion can roughly be estimated from 
longitude of periastron in 1998 and in 2007/2008 season as 
about 80 years. Further observations are, however, needed to 
precisely determine the rate of the apsidal advance and its 
quantitative interpretation in terms of the tidal distortion, 
relativistic and gravitational perturbation effects. 



7 POSSIBLE ASTROMETRIC AND DIRECT 
DETECTION 

VW LMi (HIP 54003) was astromctrically observed during 
the Hipparcos mission. No positional disturbance or accel- 
eration terms in the proper motion were found and there is 
no Double/Multiple Systems Annex flag (H59 field) in spite 
of the fact that VW LMi is member of multiple system with 
orbital period much shorter than the time span of the Hip- 
parcos astrometry. The apparent simplicity of VW LMi is 
reflected in small formal error of its annual parallax, 0.90 
mas. 

The reasons for negative detection in the Hipparcos as- 
trometry are very probably two (i) the orbital period is close 



to one year and orbital wobble mimics parallactic motion 
and changes apparent annual parallax (ii) Hipparcos astrom- 
etry refers to the photocenter of VW LMi which substan- 
tially decreases the astrometric wobble. Ecliptical latitude 
of VW LMi is about = 22.325°. The inchnation angle of 
the outer orbit is 41234 ~ 64.1±4.2° (see Section Then 
sin/3 ~ cos 41234, which means that parallactic ellipse and 
relative orbit of the photocenter would be ellipse of similar 
shape (with practically the same period of apparent mo- 
tion). The chances of disentangling the effects and solving 
the problem depend on the longitude of the ascending node. 

Outer-orbit parameters (Table[S| provide projected ma- 
jor axis as 01234 sin 11234 = 1.456±0.019 AU. Neglecting 
small eccentricity of the orbit the maximum angular sep- 
aration of the components is 7rai234. Assuming 11234 = 
64.1 ±4.2° (Section [8| the maximum angular separation is 
just about 13 mas. 

The full amplitude of the photocenter motion can be 
found assuming {Ls + L4,)/{Li + L2) = 0.42 and 034/012 = 
1.09 and a = 13 mas as: 



i^34 + iJl2 

For VW LMi we get aphot = 5.0 mas. This is compara- 
ble to the typical residuals seen in the Hipparcos astrometric 
solutions. Some improvement of the astrometric solution can 
be obtained using proper motion as known parameter from 
ground-based observations. The solving of the orbit might 
be complicated by variability induced motion (VIM) caused 
by light changes of the eclipsing pair. Therefore, VW LMi 
is rather difficult target for astrometric modelling. With ex- 
pected maximum separation of the components of about 13 
mas and with total brightness only Vmax ~ 8.08 {K = 7.20) 
it would not be easy target neither for interferometric ob- 
servations. 



8 ABSOLUTE PARAMETERS OF THE 
COMPONENTS 

Parameters of the co mponents can b e determined by a sim- 
ple procedure used in iPribuUa et aL I (|2"006 ). Using new de- 
termination of the inclination angle i = 79.0° and the pro- 
jected total mass of the contact pair (Mi -I- M2) sin'' ii2 = 
2.231 ±0.023 Mq, we obtain (Mi±M2) = 2.359 ±0.024 Mq 
0. The outer, 355-days orbit defines the mass ratio for the 
two pairs, (Mi ±M2)/(M3 ± M4) = 1.074±0.025. Therefore, 
the true (not the projected) mass of the second spectroscopic 
binary is (M3±M4) = 2.196±O.O73Af0 . Using the projected 
mass (M3 ± Af4)sin^i34 = 1.785 ± 0.011 Mq, we estimate 
the inclination of the orbit of the second pair to be about 
68.9 ± 1.5°. The outer, 355-day orbit is even less inclined to 
the sky: with (Mi±M2 + M3±M4) = 4.56±0.07 Mq and the 
projected total mass of only 3.32±0.10 Mq, we obtain 21234 
= 64.1±4.2°. In the view of large uncertainties, the planes 
of the detached pair and mutual wide orbit could still be 
coplanar. 

Individual masses of components are then Mi = 1.66 

^ we neglect error of the inclination angle depending on spectro- 
scopic determination of the third light which can systematically 
be wrong 
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Table 6. Long-term evolution of the non-eclipsing binary orbit. Parameters of contact-binary orbit and wide mutual orbit were adopted 
from the global solution in [5] Standard error of the last digit is given in parenthesis. 



Parameter 1998 1999-2000 2002 2004-5 2007-8 

HJD (2 400 000-I-) 50852-50960 51261-51673 52277-52391 53060-53836 53823-54650 

# of points 40 19 73 71 36 

634 0.032(6) 0.030(6) 0.032(5) 0.038(4) 0.039(4) 

W34 [rad] 1.43(19) 1.81(12) 2.04(15) 2.17(14) 2.18(9) 

K3 [km s-i] 63.69(37) 62.72(27) 65.00(39) 63.72(28) 63.59(34) 

K4 [km s-i] 65.67(42) 66.47(32) 65.51(43) 65.22(33) 65.79(39) 



M0, M2 = 0.70 Mq, M3 = 1.11 Mq, and A/4 = 1.09 M©. 
Masses of the components are rather inconsistent with the 
fact that all four components show practically the same spec- 
tral type. Assuming L ~ M'^'^^ and that all four components 
produce energy as MS stars results in Z/34/I/12 ~ 0.49. In the 
region of Mg I triplet we arrived at rather compatible value. 
The only remaining problem is the fact that the third light as 
found from photometry is virtually the same in the B and V 
passbands indicating the same Te/ / for both binaries which 
is incompatible with their masses. The second binary should 
have components of about GOV spectral type according to 
their masses. Due to the fact that we have spectra just in 
short wavelength region, the individual component spectra 
cannot be reliably disentangled. The issue could be resolved 
either by echelle spectroscopy and/or multi-color photome- 
try. 



9 CONCLUSIONS 

Tightest known quadruple system VW LMi is really unique. 
It definitely deserves additional observations and analysis. 
The study of this system could bring light on (i) the evolu- 
tion and origin of binary stars in multiple systems of stars 
(ii) tidal interaction of third body and its influence on the 
Roche geometry (iii) long-term evolution of orbits in tight 
multiple systems. The system is useful for further analysis 
since all four components are visible in the spectrum. The 
determination of the individual component's parameters like 
T^ff, log g, metalic i ty cou ld benefit from spectra disentan- 
gling (see 'Hadraval (|l994l )). The absolute parameters of all 
components which had to originate at practically same time 
could be used to determine the evolutionary status and his- 
tory of the system. 

The study of VW LMi is, however, complicated by (i) 
the outer orbital period being close to one year which com- 
plicates its phase coverage by Earth-bound observer (ii) very 
small angular separation of the components making direct 
resolving of the mutual wide orbit and reliable determina- 
tion of 11234 difficult (which would definitely improve de- 
termination of individual masses), (iii) fast rotation of the 
components of the contact binary making RV measurements 
unsure and rectification to real continuum impossible (due 
to the line-blanketing) (iv) eclipses and Roche geometry in 
the contact pair making usual assumption of the spectra 
disentangling techniques not valid (line profiles of individ- 
ual components cannot change with orbital revolution). 

It is also interesting to note, that VW LMi and 
HD95606 show very similar proper motion, parallax and RVs 



(see Table [21) - the stars definitely form a loosely bound pair 
and all components v ery probably evolved from the same 
protostellar cloud (see lOswalt et al.| HqOT.)). 

The orbit of the detached pair in VW LMi with 7.93- 
days period is almost circular. This is the case of another two 
quadruples detected by the DDO ob servations, TZ Boo an d 
V2610 Oph (DDO series No. XIV, iPribulla et al. 
That means that either it evolved with such orbit or it was 
circularized by the gravitational interaction with the contact 
binary. Detected apsidal motion in VW LMi requires further 
observing to reliably determine apsidal period. Investigation 
of the observed eccentricities of the second binaries and their 
predicted synchronization timescales in quadruple systems 
with contact binaries could shed li ght on the age and evolu- 
tionary status of contact binaries (|Rucinskill2008l ) . 

Better characterization of VW LMi calls for long-term 
monitoring to cover the whole 355-days orbital cycle. Es- 
pecially, times of minima should be obtained free of any 
systematic effects. The understanding of the system would 
greatly benefit from visual orbit obtained by means of 
long-baseline interferometry. Multi-color photometry and/or 
echelle spectroscopy could lead to reliable determination of 
component's temperatures and luminosities. 
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